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The linear response of synchronized time-delayed chaotic systems to small external perturbations, 
i.e., the phenomenon of chaos pass filter, is investigated for iterated maps. The distribution of 
distances, i.e., the deviations between two synchronized chaotic units due to external perturbations 
on the transfered signal, is used as a measure of the linear response. It is calculated numerically 
and, for some special cases, analytically. Depending on the model parameters this distribution has 
power law tails in the region of synchronization leading to diverging moments of distances. This 
is a consequence of multiplicative and additive noise in the corresponding linear equations due to 
chaos and external perturbations. The linear response can also be quantified by the bit error rate of 
a transmitted binary message which perturbs the synchronized system. The bit error rate is given 
by an integral over the distribution of distances and is calculated analytically and numerically. It 
displays a complex nonmonotonic behavior in the region of synchronization. For special cases the 
distribution of distances has a fractal structure leading to a devil's staircase for the bit error rate 
as a function of coupling strength. The response to small harmonic perturbations shows resonances 
related to coupling and feedback delay times. A bi-directionally coupled chain of three units can 
completely filtered out the perturbation. Thus the second moment and the bit error rate become 
zero. 

PACS numbers: 



I. INTRODUCTION 

Chaotic systems which are coupled to each other, can 
synchronize to a common chaotic trajectory [TH3] ■ Even 
for large delay times of the exchanged signals the dy- 
namic systems can completely synchronize without any 
time shift [5 10J. This phenomenon is presently attract- 
ing a lot of attention, partly because of its counter- 
intuitive fundamental aspect of nonlinear dynamics, and 
partly because of its potential for secure communication 
with chaotic signals [TTJfn]. In fact, broadband commu- 
nication with synchronized chaotic semiconductor lasers 
has recently been demonstrated over 120 km in a public 
fiber network (15) . 

Several methods for secure chaos communication over 
a public channel have been suggested [16]. One of these 
utilizes a phenomenon which was coined "chaos pass fil- 
ter" [T7] [TS] . A chaotic receiver which is driven by the 
chaotic trajectory of a sender plus a (small) message re- 
sponds essentially to the trajectory but not to the mes- 
sage. Thus the chaotic system filters out any perturba- 
tion and it is possible for the receiver to recover the mes- 
sage by subtracting its own dynamics from the incoming 
chaotic signal. 

Chaos pass filter is observed in experiments on elec- 
tronic circuits and lasers and in simulations of chaotic 
systems. But the underlying physical process is still not 
well understood. When coupled chaotic units synchro- 
nize, their chaotic trajectories are attracted towards the 
synchronization manifold, thus the synchronized system's 
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dynamics is restricted to the synchronization manifold. 
For stable chaos synchronization any random perturba- 
tion perpendicular to the manifold will exponentially de- 
cay to zero. However, this does not necessarily mean 
that a permanent perturbation like a message is damped 
as well. 

In this paper the phenomenon of chaos pass filter, i.e. 
the linear response of synchronized chaotic systems to an 
external perturbation is investigated. Research on linear 
response of general chaotic systems has been done before, 
mainly in the context of linear stochastic systems with 
multiplicative and additive noise [121110] ■ Here, we use 
iterated maps in order to obtain analytical results. In 
some respect, for example for phase synchronization |21) . 
chaotic maps have different properties than chaotic flows. 
But with respect to complete synchronization, maps and 
flows are very similar and many of the obtained results 
are also observed in numerical simulations of chaotic dif- 
ferential equations. 

This article is structured as follows: In section [H] the 
model of two coupled chaotic units is introduced and 
the linearized system's equations are derived. Section [ITT] 
defines quantities which measure the linear response to 
noisy signals, such as the second moment and the bit er- 
ror rate. The previously introduced model is studied ex- 
tensively. Large excursions off the synchronization mani- 
fold are observed which lead to a continuum of diverging 
moments. Bit error rates, devil's stair cases, resonances 
and Lyapunov spectra are calculated. In section |IV| the 
investigation is extended to more complicated models, 
such as a chain of three units and a four units network. 
Finally, the last section summarizes and discusses our 
results. 
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FIG. 1: Setup of two coupled chaotic units with either uni- 
or bi-directional (dashed line) coupling. A perturbation m is 
added to the exchanged signal at unit 1, i.e., the sender. The 
transmitted signal has a time delay r. 



II. THE TWO UNITS SETUP 

The simplest model for a coupled interacting chaotic 
system comprises of just two chaotic units which can ei- 
ther be uni- or bi-directionally coupled by a function of 
their internal variables. The exchanged signal usually 
has some time delay r and in order for the bi-directional 
setup to be able to synchronize we in general include 
some self-feedback having the same time delay. The ex- 
ternal perturbation m is added to the transmitted signal 
at the sender as depicted in Fig. [I] The unidirectional 
configuration has been realized in a communication net- 
work over 120 km |15j . In terms of chaos communication, 
where such a setup is called chaos masking |16j . one can 
think of the perturbation being a message which shall be 
secretly transmitted on top of a chaotic carrier signal. 
The message is small compared to the carrier signal and 
for our purpose its content may be considered as random 
noise. Hence we are, in the following, interested in the 
linear response of the receiver to noise. 

Note that for chaos communication a uni-directional 
setup is not secure since an eavesdropper who, for ex- 
ample, knows all the details can use an identical copy 
of the receiver, synchronize it as well and extract the 
message. In contrast, two chaotic units which interact 
by bi-directional transmission have an advantage over an 
attacker driven by a uni-directional signal [22] - 

The dynamics of these two chaotic systems is given by 
the following sets of iterated equations 

• uni-directional setup 

x t +i = (1 - e)/(a? t ) + e/(a; i _ T ) 
Vt+i = (i - e)/(y*) 

+ eKf(yt-r) + e(l - K)f(x t - T + m t _ T ) , (1) 



• bi-directional setup 

xt+i = (1 - e)f(x t ) 

+ enf{x t - T ) + e(l - K)f(y t - T ) 
Vt+i = (1 - e)/(y t ) 

+ enf(y t - T ) + e(l - K)f(x t - T + m t _ T ) , 



(2) 



where t is a discrete time step. The parameter e measures 
the contribution of all delay terms while parameter k 
measures the relative strength of the self-feedback. They 
are chosen such that the dynamics, Xt and yt, always 
stays in the unit interval [0, 1]. For the function f(x) one 
of the following discrete maps is used: 
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FIG. 2: Phase diagram for two coupled Bernoulli units with 
a — 1.5. System is completely synchronized in regions II + 
III (- -) for the the uni-directional setup and in regions I + II 
(-) for the bi-directional setup, see Eq. |6]| and Eq. Q. 



Bernoulli map 

f(x) = (ax) mod 1 



• Tent map 

/(*) = 

• Logistic map 



i x for < x < a 

^ (1 - x) for a < x < 1 



f{x) — a x(l — x) 



(3) 



(4) 



(5) 



Without noise, m t = 0, the spectra of Lyapunov ex- 
ponents and hence the phase diagram of synchronization 
can be calculated analytically for the Bernoulli system in 
the limit of large delay, r —> oo [52 [H] . The parameter 
range for which the system in the stationary state syn- 
chronizes completely, i.e., x t — yt, is depicted in Fig. [2] 
and is given by the inequality 



• uni-directional setup 

l-a(l- e) 



k < 

• bi-directional setup 
a-1 



ac 



2ae 



< k < 



1 + 2ae 



2ae 



(G) 



(7) 



The effect of small noise can be calculated by lineariz- 
ing equations ([TJ) and |2| in the vicinity of the synchro- 
nization manifold x t — yt = s t . The noise leads to a 
small deviation d t — yt — x t which is determined by the 
linear equations 



uni-directional setup 

dt+i = (1 - e)/ t 'rfi + enft_ T d t - 



(8) 
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• bi-directional setup 
d 



H+l 



= (1 - e)f t d t + e K fl_ T d t _ T 

+ e(l - K)f t _ T (m t - T - dt-r) 
= (1 - e)f t d t + (2k - l)e/ t '_ T d t _ T 
+ e(l ~ K)f't-T m t-r 



(9) 



where // is the derivative of f(x) at the synchronized 
trajectory St- These equations determine the region of 
synchronization. Without noise dt decays to zero for 
stable synchronization and increases exponentially oth- 
erwise. Note that in this case, i.e., without any noise 
present, these equations can be solved analytically for 
Bernoulli units, where / t ' = a = const, and determine 
the phase diagram of Fig. [2] 



III. LINEAR RESPONSE 

We want to understand how a synchronized chaotic 
system responds to small perturbations, i.e., the linear 
response of the system. In the following we analyze the 
linear response by means of the moments of the distri- 
bution of dt , the bit error rate of a transmitted random 
message and resonances of harmonic perturbations. 



A. Moments 

Without an external perturbation the coupled system 
can synchronize completely, see Fig. [2] With noise how- 
ever the deviation from the synchronization manifold d t 
does not decay to zero but has a distribution around zero. 
We are interested in the stationary distribution which 
develops after a transient time and which is caused by 
the transverse Lyapunov exponents. This distribution is 
characterized by its moments which are defined as 



Xn 



lim 



<M)^o (H«) ' 



(10) 



where (. . .) means an average over the distribution of 
m and d, respectively and n denotes the order of the 
moment. 

For investigating the moments we model the external 
perturbation m t as random numbers with a uniform dis- 
tribution in the interval [— M, M\. 

In order to obtain analytical results let us consider the 
simplest case, the uni-directional setup without any time 
delay, r = 0. Using the substitution a = e(l — k) the 
system reduces to a simple master-slave setup without 
self-feedback and the linearized equation ^ takes the 
simple form 



H+l 



(1 ~ a )ft d t + af' t m t 



(11) 



For Bernoulli maps f' t is a constant and for the other 
chaotic maps we assume f t to be an uncorrelated random 
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FIG. 3: Distribution of jdj for two uni-directionally coupled 
logistic maps with a = 4, k = 0, r = and M = 10 -8 in a 
log-log plot. Dashed red line shows a power law fit. There 
is a cut-off from the power law behavior for small d due to 
the additive noise which prevents the system from perfectly 
synchronizing. 



number. In the latter case Eq. (11) is a discrete linear 



stochastic equation with multiplicative and additive noise 
of the form 



Xt+i = JtXt + Vt 



(12) 



In the theory of discrete linear stochastic equation it is 
well known that such an equation may lead to stationary 
distributions which have a power law of the form [5D1 



p{x) 



1 



where fi is a solution of the following equation 

(y- 1 ) = i. 



(13) 



(14) 



As simulations show the distribution p[d) may indeed 
follow a power law. If the factor (1 — a)|/ t '| in Eq. (11) 



takes on values larger than one, dt can temporarily ex- 
plode, i.e., the system can have very large excursions from 
the synchronized state, leading to a power law. For the 
logistic map with a = 4, for which the distribution of /' 
is given by 



P(f) = 



1 



(15) 



the maximum slope is |/'| = 4 and hence we expect to 
find large excursions of d ( for e < | . Fig. [3] shows a 
typical distribution of \d\ for two uni-directionally cou- 
pled logistic maps. The distribution has a peak at values 
slightly larger than the maximum noise strength M, since 
the additive noise in the systems prevents it from being 
perfectly synchronized. For larger distances the distribu- 
tion follows, to a good approximation, a power law. 

We did not succeed in calculating the distribution 
p(d) analytically, but the linear response can be derived. 
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FIG. 4: Second moment \2 (blue points) and cross correlation 
C (red squares) for two uni-directionally coupled logistic maps 
with a — 4, r = 0, M = 10 -8 and k — (such that a = e) as 
a function of e. Dashed curve shows analytical results for \2 
whereas other results were obtained from simulations. 



Squaring Eq. (Illl) one finds for the second moment 



« 2 (/' 2 > 



(d 2 ) 

(to 2 ) l-(l-a) 2 (/' 2 ) 



(16) 



The second moment diverges when the denominator be- 
comes zero meaning that the distribution p(d) has devel- 
oped a power law tail or in other words that rare but 
large excursions from the synchronization manifold oc- 
cur. For the Bernoulli map with a = 1.5 the second 
moment diverges for a < a m — | which coincide with 
the synchronization transition given by 



ln(l-a) + (ln|/'|) <0. 



(17) 



Without multiplicative noise all moments diverge at 
the same parameters. For the logistic and the tent map, 
however, the second moment diverges already inside the 
region of synchronization. For the logistic map with a = 
4 and the mean values (/"•) = 8 and (ln|/'|) = In 2 
which can be analytically computed from the distribution 
Eq. (151 we find the synchronization transition to be at 
a s = "j and the second moment to diverge at a m = 1 — 
0.646. For the tent map with a = 0.4 where (/' 2 ) = 
(ln|/'|) re 0.67 we find a s 0.49 and a m = 

0.51. Fig. j^J shows the second moment %2 
together with the cross-correlation C, which is a measure 
for the synchronization, for a system of coupled logistic 
maps. The numerically obtained results for the second 
moment agree with the analytical results Eq. ( 16 ) and 
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Vs 
^ and 

1 - 



one can see clearly that the second moment diverges even 
though the system is still synchronized. 

Fig. [5] shows a trajectory of d t for logistic maps where 
the coupling parameters are once chosen such that the 
second moment exists and once that it diverges. In both 
cases the system is synchronized, i.e., C = 1, but in case 
of a diverging second moment large excursions from the 
synchronization manifold occur. On the contrary, for a 
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FIG. 5: Trajectory of dt for two uni-directionally coupled 
Logistic maps with a = 4, k = 0, r = 0, M = 10 -8 and (a) 
e = 0.6 and (b) e = 0.7. In both cases the cross-correlation is 
C = 1 but for (a) large excursions from the synchronization 
manifold occur and the second moment diverges, \2 ~ 4632, 
whereas for (b) the second moment is finite, x% — 14. Inset 
in (b) shows that the deviations from the synchronization 
manifold is in the order of the magnitude of the noise. 



finite second moment the deviations from synchroniza- 
tion are mainly determined by the additive noise term in 
Eq. ( 11 ) and are of the magnitude of M. 



Note that we obtain similar plots to Fig.[3][5]for coupled 
tent maps but due to the limited space we restrict the 
presentation of our results at this point to coupled logistic 
maps. 

Power law tails of the distribution of a stochastic pro- 
cess with multiplicative and additive noise have been dis- 
cussed in the context of chaos synchronization, before [2] • 
This phenomenon has been called "on-off intermittency" 
and was discussed in the vicinity of the synchronization 
transition. Here the phenomenon is triggered by the mes- 
sage transmitted via the chaotic signal, and it is observed 
deep inside the region of synchronization. 



The linear equation (11) also determines higher mo- 
ments of p(d). By taking the nth power of Eq. (11) we 
find that \n diverges at a coupling e n given by 

1 = (1 -e n ) n (/' n ). (18) 

Hence, for Bernoulli maps all moments diverge at the 
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FIG. 6: The threshold e n below which \ n diverges for two 
urii-directionall y co upled logistic maps with a = 4, n — and 
t = 0, see Eq. (pi. 



synchronization threshold e s . For logistic maps, however, 
all moments diverge at a different coupling strengths e n 
given by 



2T( , - 



(19) 



Starting from a strongly synchronized state and decreas- 
ing e the distribution p(d) broadens and eventually fol- 
lows a power law. The different moments of the distribu- 
tion successively start to diverge starting with the highest 
one Xoo and continuing with lower and lower ones until fi- 
nally for n — ¥ the synchronization threshold is reached, 
where (ln|d|) diverges. Figure [6] shows e„ as a function 
of n. Note that the broadening of the distribution by de- 
creasing e is related to a decreasing power law exponent 
p. Comparing equation (18) and (14) we see that p is 



related to the order of the moment n by p, = n + 1. 

Up to now we have discussed the uni-directional setup 
only. The linearized equations for the bi-directionally 
coupled system without time delay are very similar and 
with the substitution a = e(l — k) read 



d t+1 = (1 - 2a)f t d t + af' t m t 



(20) 



However this system is harder to analyze analytically 
since due to the mutual interaction the distribution of 
/' is rendered, p(f') does not correspond to the distri- 
bution of an isolated unit anymore as it was the case for 
the master-slave setup. Only for Bernoulli maps where 
/' = const the distribution p(f') is not affected by the 
mutual coupling. 

The second moment is given by 



(d 2 



2 (/' : 



(m 2 ) l-(l-2a) 2 (/' 2 ) 



(21) 



and it diverges if the denominator approaches zero. Al- 
though the mutual interaction renders the distribution of 
/' and therefor we cannot compute the mean (/' 2 ) nor 
(In | /'I) for the tent and logistic map, respectively, we 
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FIG. 7: Second moment \2 (blue points) and cross correlation 
C (red squares) for two bi-directionally coupled (a) Bernoulli 
and (b) logistic maps with a = 1.5 and a = 4, respectively, 
as a function of e. Other parameters are r = 0, M — 10~ 8 
and ft = (such that a — e). Dashed curve shows analyt- 
ical results for X2 whereas other results were obtained from 
simulations. 



obtain a surprisingly good agreement with numerical re- 
sults if we use, as a first approximation, the distribution 
of a single unit for computing the means. Fig.[7]shows \i 
obtained from numerical simulations together with ana- 
lytical results for Bernoulli and logistic maps. Only at 
the transition from finite to diverging second moments 
the simulations deviate slightly from the analytical re- 
sults in case of logistic maps. This indicates that p(f') is 
changed only moderately by the bi-directional coupling. 
Closer investigations of p(f') have indeed confirmed that 
the distribution is hardly changed for r = 0, see for ex- 
ample Fig. [8] (a) which shows p(f') for the receiver of 
two synchronized bi-directionally coupled logistic maps 
together with the distribution of a single logistic map 
given by eq. ( 15 ). 



Note that, as for the uni-directional setup, the second 
moment already diverge within the region of synchro- 
nization for logistic and tent maps whereas for Bernoulli 
maps the second moment diverges at the synchroniza- 
tion transition. With increasing coupling the effect of 
the noise becomes stronger and, hence, the susceptibility 
is not symmetric around a — 0.5. 

In the general case described by equation ^ and ^ 
for the uni and bi-directional setup, respectively, a time 
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FIG. 8: Probability distribution p(f') of the receiver of 
coupled synchronized logistic maps for (a) undelayed bi- 
directional setup with e = 0.4 and r = and (b) time-delayed 
uni-directional setup with e = 0.7 and r = 100. Other param- 
eters are k = and M = 10 -8 . Dashed (red) line sho ws t he 
distribution of a single logistic map given by equation ( 15 I. 
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FIG. 9: Auto correlation A(t) of dt for two uni-directionally 
coupled logistic maps with a — 4, e = 0.7, K = 0, r = 1000 
and M = 10~ 8 . 



delay r is present. The trajectory d t shows distinct auto- 
correlations at integer multiples of r see Fig. [9] Un- 
fortunately, we generally do not know the magnitude of 
the correlation between dt and rft_ T , which also changes 
for varying coupling parameters. Additionally the distri- 
bution of /' is altered tremendously for a time delayed 
system compared to the distribution of a single unit, see 
Fig.[8](b). Hence we cannot calculate the second moment 
analytically but have to rely on numerical simulations. 

A numerical investigation of the time delayed system 
shows that the behavior of the moments is similar to a 
system with r = 0. For Bernoulli maps all moments 
diverge at the synchronization threshold whereas for lo- 
gistic and tent maps \2 already diverges inside the re- 
gion of synchronization. Numerical results for two uni- 
directionally coupled logistic maps are exemplary shown 
in Fig. [lO] as a binary plot where the system is assumed to 
be synchronized when the cross correlation exceeds some 
threshold 9 S and the second moment is assumed to be 
finite when it is smaller than some threshold 9 m . 



FIG. 10: Phase diagram for two uni-directionally coupled lo- 
gistic maps with a — 4, r = 100 and M — 10 -8 . Gray regime 
shows synchronization region, i.e., area where C is larger than 
a threshold 8 C — 0.999. Black regime shows region where X2 
is finite (and the system is synchronized), i.e., X2 smaller than 
a threshold 9 m — 50. 



B. Bit Error Rate 

The bit error rate (BER) measures the quality of the 
transmission of a (binary) message. For a proper recon- 
struction of the transmitted message it is crucial that 
the system is synchronized. However, the message is an 
external perturbation to the system that potentially de- 
stroys the synchronization. Hence the BER tells us how 
an external perturbation influences the synchronization, 
i.e., it is an indirect measure of the linear response. 

In the following the sender transmits a binary message 
of the form m t = ±M with Af <1 and (m) — which is 
reconstructed by the receiver by subtracting its own state 
from the received signal |16) . The recovered message rht 
is then given by 



rht = i x t + mt) - y t = m t - d t . 



(22) 



If both units were perfectly synchronized, i.e., dt = 0, 
the original message would be perfectly recovered by the 
receiver. For a successful reconstruction of a binary mes- 
sage it is sufficient that both messages, mt and rht, have 
the same sign, m t fht > 0. 

The quality of the reconstruction is given by the BER 
r which is defined as the the number of incorrectly re- 
covered bits normalized by the total number of received 
bits, 



# incorrectly recovered bits 
# received bits 



(23) 



Note that by just guessing the bits we would be correct 
in half of the cases on average and hence obtain r = 0.5. 

The BER is related to the distance d by an integral 
over the distribution p{d). If the absolute value of d is less 
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than the absolute value of the messages, \d t \ < M, the 
message is recovered correctly with probability 1 because 
the sign of the message is not changed. In the other 
cases, \dt\ > M, the sign of the message is changed with 
probability |, thus the BER is given by 



(24) 



r= \ I J P (d)Ad + J p{d)dd 

\— 00 M 

= i (l- J P (d)ddj . 



Note that this definition of the BER assumes that the 
system has relaxed to a stationary distribution p(d). In 
fact, this definition may be only an upper bound since 
lower BERs may be achievable if one uses additional in- 
formation about the transmitted signals |27J. 

The distribution p(d) is known analytically only in 
some special cases, see Appendix. In general, the BER 
has to be determined by means of computer simulations. 
Figure [TT] shows results from simulations for the uni- and 
bi-directional setup without time delay for Bernoulli, tent 
and logistic maps. For an unsynchronized system the 
BER is at its maximum r = 0.5. In contrast to the sec- 
ond moment which becomes finite at the synchronization 
transition only for Bernoulli maps, the BER drops down 
to smaller values for all maps as soon as the system syn- 
chronizes. That means that the BER is smaller than 
r = 0.5 although d has large excursions from the syn- 
chronization manifold and the second moment diverges. 
Hence the sign of the deviation has correlations to the 
original message even though its amplitude has still a 
broad distribution. In this parameter range information 
can already be successfully transmitted. Although the 
BER is very high, one may apply methods from informa- 
tion theory to derive the message. 

The BER shows a staircase structure for Bernoulli and 
also for tent maps. As derived in the Appendix this 
devil's staircase is related to a fractal structure of the 
distribution p(d). For Bernoulli maps we find a broad dis- 
tribution of d for small e, see Fig. 12 'a). With increasing 



e the distribution obtains more and more a fractal struc- 
ture while the support stays connected, see Fig. [l2^b), 
until eventually the distribution changes to a peaked 
structure with a fractal support, see Fig. [l2](c). For the 
uni-directional setup the peaked distribution occurs for 



e > | and the BER locks into rational values r 



with k and q being natural numbers and ~ < 
For the bi-directional setup the distribution is peaked for 
| < e < I and the BER is zero for | < e < | since the 
absolute values of the distances are less than the message 
amplitude, \d\ < M, in this interval. Note that the frac- 
tal properties of the distribution p{d) are related to the 
theory of iterated function systems [28]. For Bernoulli 
maps Eq. (Ill and ( 20 1 give iterations of two linear func- 
tions with m t = ±1. Iterating a few randomly chosen 
functions can lead to fractal distribution. For details see 



Appendix. 

The distribution p(d) for coupled tent maps has a sim- 
ilar behavior to the one for Bernoulli maps. It also shows 
a peaked structure which is related to a staircase in the 
BER, see Fig. 12 d) . But, in contrast to Bernoulli maps, 
the distribution for the tent map system can have very 
long tails, see Fig. [l2^e). These occur in the parameter 
range where the system is synchronized but the second 
moment diverges. For Bernoulli maps there is no such 
parameter range since both transitions take place at the 
same point. 

For logistic maps the distribution p(d) can also have 
power law tails. But unlike to Bernoulli and tent maps it 
does not show a peaked structure. The distribution has 
always a connected support due to the broad distribution 



of the multiplicative noise /', see Fig. 12 f). 



The BER for a system with (large) time delay is ex- 
emplary shown for bi-directionally coupled tent maps in 
Figure |11| The synchronization transition is also indi- 
cated in the plot and illustrates that the BER, similar 
as for the undelayed system, is reduced as soon as the 
system synchronizes. 

Note that our results show that the integral over the 
distribution of distances p{d) , Eq. p4| , which determines 
the BER, is insensitive to long tails of p(d). Although the 
responses Xn diverge at different coupling strengths e n , 
the BER is smooth as a function of e. 

In our model the BER is rather high compared to re- 
ported values of analog and other digital systems. Hence 
our investigation just gives a qualitative explanation of 
chaos pass filter. But note, that we are using bits of 
length L = 1 for our investigation. The bit error can be 
exponentially reduced by increasing L. 



C. Resonances 

In the previous subsections we investigated the lin- 
ear response of synchronized chaotic units to a random 
perturbation mt- We in particular analyzed the dis- 
tance d of the units given by the linearized equations 
^ and Q for the uni and bi-directional setup, respec- 
tively. For Bernoulli maps these linear equations have 
constant coefficients and hence any arbitrary perturba- 
tion can be decomposed into its Fourier modes. In this 
special case it is sufficient to investigate the linear re- 
sponse of the system to a harmonic perturbation of the 
form mt = Mexp(— iut). The system responds with the 
identical frequency d t = Dexp(— iut), and the amplitude 
of the recovered signal m = mt — d t is given by 



\M - D\ , 



(25) 



with the complex amplitude D. We obtain the following 
results 



• uni-directional case 



D = M- 



e(l — n)a 



- (1 - e) 



ae 



■ ena 



(26) 




FIG. 11: Bit error rate r for uni- and bi-directionally coupled (a) Bernoulli (b) tent and (c) logistic maps with a = 1.5, a = 0.4 
and a = 4, respectively. Other parameters are r = 0, k = and M = 10 -8 . The synchronization transition is indicated by e a 
and the transition to a finite second moment by e m . 
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(e)Tent map, bi-directional, e = 0.25 (f)Logistic map, bi-directional, e = 0.1 

FIG. 12: Probability distribution p(d) for different setups with parameters r = 0, k = and M = 10 -8 . 



bi-directional case 



D = M- 



e(l — k)<2 



j(t+i) - (1 - e)ae- JWT - e(2« - l)a 



(27) 



Figure [M] shows |m| as a function of w close to the 
phase boundary of synchronization. The amplitude of the 
reconstructed signal shows peaks (resonances) separated 
by a distance 27r/ (t+1) which are caused by the first term 
of the denominator in equation ( 26 ) . The magnitude 



of the peaks is modulated due to the discrete nature of 
the systems equations. The resonances diverge at the 
synchronization boundary. 

For the logistic and tent map the coefficients of the 
linear equations ^ and ^ depend on time. Thus an 
exact Fourier decomposition of the transmitted signal is 



not possible. However, our numerical results show that 
the corresponding Fourier component of the response frit 
shows resonances as well. The power spectrum of the 
recovered message has a clear peak at the frequency of 
the harmonic perturbation, see Fig. |15| Surprisingly, no 
higher harmonics are observed. 

These results show that a chaotic system can func- 
tion as a sharp harmonic filter. The transmitted signal 
is irregular and the harmonic perturbation is arbitrarily 
small. Nevertheless, the receiver can filter out this per- 
turbation with high precision. It would be interesting 
to investigate this harmonic filter with respect to secret 
communication. In fact, the response of synchronized 
chaotic semiconductor lasers to a harmonic perturbation 
has been investigated in Ref. [25] . 
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FIG. 13: Bit error rate r for a bi-directional tent map system 
with delay r = 100. The black solid line indicates the syn- 
chronization transition, i.e., the boundary for which all cross 
correlations are larger than a threshold, C > 9 S — 0.999. The 
other parameters are a = 0.4 and M = 10~ 8 . 




FIG. 14: Amplitude of the reconstructed signal m for different 
signal frequencies cj. Bernoulli map, uni-directional coupling 
with e = 0.8, k = 0.55 and r = 50. 



D. Transverse Lyapunov Spectra 

For stable chaos synchronization the distance dt decays 
to zero such that the system's dynamics is restricted to 
the synchronization manifold. An external perturbation, 
however, drives the system away from this manifold. The 
competition between these two mechanisms results in the 
linear response investigated in the previous sections. 

The relaxation to the synchronization manifold is de- 
scribed in terms of transverse Lyapunov exponents. For 
stable synchronization all transverse exponents are nega- 
tive on average whereas in the unsynchronizcd case pos- 
itive exponents exist. Close to the transition the max- 
imum transverse exponent becomes very small until it 
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FIG. 15: Power spectrum S(to) of a reconstructed message m 
in a semilogarithmic plot. Original message is of form mt = 
Msina>o£ with M = 10 -8 and uj = 5 ■ 2-r/t. Tent map, uni- 
directional coupling with e = 0.6, K = 0.4 and r = 100. Note 
that the underlying periodic structure is due to the sampling 
frequency. 



eventually crosses zero at the transition. This results in 
a slowing down of the relaxation to synchronization and 
in a divergence of the response close to the transition. 
Diverging moments and power law tails are in particu- 
lar related to local positive transverse Lyapunov expo- 
nents. Although the transverse exponent is negative on 
average in the synchronized case there can occur positive 
local ones which lead to temporary large excursions away 
from the synchronization manifold. These positive local 
transverse exponents may be caused by unstable peri- 
odic orbits or invariant subjects of the chaotic manifold 

G2BEI]. 

Following these arguments, one might expect that the 
qualitative behavior of the response functions x„ or the 
bit error rate r is related to the behavior of the largest 
transverse Lyapunov exponent X max . However, this is 
not true. For example, consider the simple driven sys- 



tem without delay, Eq. (11). The transverse Lyapunov 
exponent is given by (17), 



ln(l 



Mf'\), 



(28) 



and decreases monotonically with increasing coupling 
strength e. However, the bit error rate r first decreases 
with e to a minimum value before it increases again, see 
Fig- [TTJ 

For a Bernoulli system with delay r the spectra of Lya- 
punov exponents can be calculated by solving polyno- 
mial equations of order r |32j . Fig. 16 shows the trans- 



verse Lyapunov spectrum as a function of e for the self- 
feedback strength k = 0.2. Again, the qualitative behav- 
ior of the largest exponent is not related to the one of 
the bit error rate r. Except for the synchronization tran- 
sition where the maximum exponent becomes negative 
and the BER becomes less than its maximum of |, the 
qualitative behavior of the Lyapunov exponents and the 
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FIG. 16: Spectrum of transverse Lyapunov exponents A and 
bit error rate r for Bernoulli map, uni-directional coupling, 
K = 0.2, r = 50. 



BER is not related. 



IV. SMALL NETWORKS 



FIG. 17: Setup of three coupled chaotic units with either uni- 
or bi-directional (dashed line) coupling. A perturbation m is 
added to the exchanged signal at unit 1, i.e., the sender. The 
transmitted signal has a time delay r. 

• uni-directional setup 

x 1 t+1 = {l-e)f{x 1 t )+ef{x\_ T ) 
x 2 +1 = (l-e)f(x 2 )+eKf(x 2 _ T ) 

+ e(l - «) /04-t + m t-r) 

x 3 +1 = (1 - e)f(x 3 ) + e K f(x 3 _ T ) + e(l - n)f(x 2 _ T ) 

(29) 

• bi-directional setup 

x) +1 = (1 - e)f(x}) + enf{x\_ T ) + e(l - n)f(x 2 _ T ) 
x 2 +1 = (l-e)f(x 2 ) + eKf(x 2 t _ T ) 

+ 6(1 - k) (1/2 f{x\_ T + mt-r) + 1/2 f(x 3 _ T )) 
x 3 +1 = (1 - e)f(x 3 ) + eKf{x 3 _ T ) + e(l - K)f{x 2 _ T ) 

(30) 



Having discussed the linear response of a coupled two 
units system in the previous chapters, we are investigat- 
ing the linear response of more complicated setups such 
as a chain of three units and a network of four units 
in this chapter. The motivation again stems from chaos 
communication. Hence, in all setups the external pertur- 
bation is added to the signal of the " sender" and can be 
thought of as a secret message. 

The systems in this section are generally to compli- 
cated for an analytical discussion of the moments \ n or 
the BER r. Hence we mainly rely on numerical simula- 
tions in the following. 

We define the second moment for different units by Xij 
where the index stands for the combination of unit i and 
j, i.e., the distance dij — x\ — x\ is used when evaluating 
Eqn. (10). Similar, for the BER the signals of units i 



and j are compared. 



A. Chain of Three Units 

In the following the system comprises of three chaotic 
units which are arranged in a line such that the first unit 
is coupled to the second and the second to the third. 
The coupling can either be uni- or bi-directional and has 
a time delay. Optionally the units are subjected to a self- 
feedback with the same delay time. The setup is depicted 



with the coupling parameters e and k as before. 

For Bernoulli units without noise, m t — 0, the syn- 
chronization properties can be analyzed analytically by 
a master stability function method [53J [33] . In the limit 
of large delays one finds that for the uni-directional setup 
only complete synchronization can occur whereas for the 
bi-directional setup complete and sub-lattice synchro- 
nization exists [34L For the uni-directional setup the 
system synchronizes completely in the parameter regions 
II + III of Fig. [2] which is identical to the synchroniza- 
tion region of two uni-directionally coupled units. The 
bi-directional system synchronizes completely in region 
II whereas in region III only unit 1 and unit 3 are syn- 
chronized. 

In case of uni-directional coupling unit 1 imposes its 
chaotic behavior onto unit 2 which in turn imposes its be- 
havior onto unit 3. All three units become synchronized 
for the same coupling parameters, however in the unsyn- 
chronized case the cross correlations are highest between 
unit 1 and 2. For a Bernoulli system, numerical results 
for the second moment and the BER together with the 
respective cross correlation for uniformly distributed and 
binary noise are shown in Fig. 18 In contrast to the 



in Fig. 17 and the system's equations are 



two units setup, where the second moment for Bernoulli 
maps diverge at the synchronization transition, the sec- 
ond moment X13 diverges already inside the synchronized 
region. As we expect, the response is smallest for the 
combination of unit 1 and 2. The BER decrease for both 
combinations, r 12 and r 13 , at the synchronization tran- 
sition, similar to the two units setup, and is in general 
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FIG. 18: Panel (a) shows second moment Xij (blue curves) FIG. 
and cross correlation Cy (red curves) for uniformly dis- and 



19: Panel (a) shows second moment Xij (blue curves) 
cross correlation dj (red curves) for uniformly dis- 



tributed random noise, panel (b) shows bit error rate ry (blue tributed random noise, panel (b) shows bit error rate Tij (blue 



curves) and cross correlation dj (red curves) for binary ran- 
dom noise as a function of e. Solid (dashed) line shows the re- 
spective results for the combination of units ij = 12 (ij=13). 
Chain of uni-directionally coupled Bernoulli maps with pa- 
rameters a = 1.5, k = 0, r = 100 and M = 10~ 8 . 



curves) and cross correlation dj (red curves) for binary ran- 
dom noise as a function of e. Solid (dashed) line shows the re- 
spective results for the combination of units ij = 12 (ij =13,). 
Chain of bi-directionally coupled tent maps with parameters 
a = 0.4, k = 0.3, t = 100 and M = 10" 8 . 



lowest for r±2- A devil's staircase occurs as well at which 
rj.3 becomes the lower BER. 

For the bi-directional setup unit 1 and 3 synchronizes 
perfectly, hence the distance di 3 and also the second mo- 
ment X13 is zero. A perturbation of the signal which 
unit 1 sends to unit 2 is completely filtered out and does 
not affect unit 3. This can easily be seen in the system's 
equations (30). Unit 1 and 3 obtain the same external in- 



put from unit 2 and since the units are identical with the 
same internal dynamics they synchronize perfectly. Note 
that unit 2 does not necessarily need to be synchronized 
in order for perfect synchronization between unit 1 and 
3 to occur. It acts as a relay which transmits the signals 
between the two units such that they can synchronize. 
Comparing the outgoing signal of unit 1, i.e. its internal 
dynamics plus the message m, with the outgoing signal 
of unit 3, the message can be recovered perfectly without 
any errors, thus the BER is r%3 = 0. 

Fig. 19 shows the second moment and the BER to- 



gether with the cross correlations for a system of tent 
maps as a function of the coupling e. The coupling con- 
stant k was chosen such that sub-lattice synchronization 



occurs for some values of e. Sub-lattice synchronization 
exists in the region where C13 = 1 but C12, C23 < 1. At 
the transition to sub-lattice synchronization X13 and ri3 
becomes zero. At the transition to complete synchroniza- 
tion the BER r 12 decreases whereas the second moment 
X12 is still diverging close to the boundary. But neither 
the BER 7*12 nor the second moment X12 becomes zero. 
The external noise prevents the trajectories, x\ = x$ and 
X2, of perfectly synchronizing. It causes the trajectories 
to deviate by a factor of the noise strength. 

In terms of chaos communication, unit 2 which obtains 
the signal of unit 1 and 3 can recover a message, which is 
added on top of the transmitted signal of unit 1, without 
any errors. In the same fashion more units can be added 
to the system in a star like setup with unit 2 being the 
relay. The additional units can also add a secret message 
on top of their transfered signals. Fig. [20| shows the setup 
for a system with four units, but in principle an arbitrary 
number of units can be added. All units of the star, apart 
from unit 2, synchronize perfectly since they receive an 
identical driving signal. Unit 2 compares the signal from 
unit 3 with the incoming signals of units 1 and 4 and 
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FIG. 20: Setup of a star like configuration with unit 2 in 
the middle. Unit 2 obtains the unperturbed signal of unit 
3, and the signal of unit 1 and 4 where a message m and 
w, respectively, is added. Unit 1, 3 and 4 can synchronize 
perfectly. Hence unit 2 can obtain m and w by subtracting 
the (unperturbed) signal of unit 3 from the signal of unit 1 
and 4, respectively. 



perfectly recovers both secret messages. Thus the hub of 
the star can simultaneously decode any number of secret 
messages. 



B. Four Units Network 

In the following we investigate a ring of four coupled 
units. For Bernoulli units the synchronization properties 
for the unperturbed system can be calculated analyti- 
cally. The stability of synchronization is determined by 
the eigenvalues of the adjacency matrix which describes 
the coupling of the units. For an adjacency matrix with 
the row sum equals one, meaning that all incoming sig- 
nals are normalized, and in the limit of large delays one 
finds that a spectral gap between the largest eigenvalue 
of 71 = 1 and the second largest eigenvalue 72 is cru- 
cial for the stability. Complete synchronization is only 
possible in the limit of weak chaos if the spectral gap 
is nonzero [24 . For a simple ring network without any 
self-feedback, where all eigenvalues are 7 = 1, no eigen- 
value gap exists and it cannot synchronize. Adding an 
additional link with the coupling strength a changes the 
eigenvalues such that a gap occurs so that the system is 
able to synchronize |35j . 

The perturbed system, where a noise m is added onto 
all outgoing signals of unit 1, is depicted in Fig. 31 and 
is described by following equations 



x t+l 


= (1- 






-r) 




x t+l 


= (1- 




- tf( x l 




- rnt-r) 


T 3 

X t + 1 


= (1- 












+ e( 




+ mt-T 


) + 


(l-a)f(x 2 t _ r )) 


X t+1 


= (1- 


e)/(4) "t 


- e/(a£ 


-r) 





(31) 



For such a network the eigenvalue gap in the unperturbed 
case, and therefor the synchronization ability, is maximal 
for a w 5/8 = 0.625. Note that for weak chaos and hence 
for complete synchronization in this network, the chaos of 
the single units has to be small. For Bernoulli maps one 
finds for the critical coupling at which synchronization 




FIG. 21: Ring network of four units plus an additional link 
of strength a in order for the network to synchronize. A 
perturbation m is added to the outgoing signals of unit 1. 



occurs 



Cr > 



1 - 1/a 
I-I72I 



(32) 



with a the parameter of the Bernoulli map. Thus for 
the maximum eigenvalue gap with a = 0.625 the system 
synchronizes for a < 1.16. 

For a Bernoulli system the second moment Xij an d the 
BER rij for combination of unit 1 with unit j = 2, 3, 4 
are shown in Fig. 22 together with the respective cross- 
correlations. Surprisingly the combination of unit 1 and 
2 has the highest second moment and BER whereas unit 
4, which is only indirectly driven by unit 1, has the lowest 
second moment and BER. 



V. SUMMARY 

The linear response of a time-delayed chaotic system 
to small external perturbations has been studied. The 
external perturbation drives the system away from the 
synchronization manifold whereas the dynamics of the 
system, quantified by transverse Lyapunov exponents, 
relaxes the system back to the manifold. The compe- 
tition between these two mechanism results in the linear 
response studied in this work. 

This investigation is motivated by chaos communica- 
tion, where a secret message is added on top of an ex- 
changed signal between synchronized chaotic units, thus 
perturbing the system. Nevertheless, the receiver is able 
to recover the secret message by subtracting its own 
chaotic trajectory from the received signal. This mecha- 
nism has been named chaos pass filter since the receiver 
filters out any external perturbation and essentially re- 
sponds to the unperturbed chaotic trajectory. 

However, our numerical and analytical investigations 
of iterated maps show that the mechanism of chaos pass 
filter is much more complex. Perturbations are not just 
damped, instead the response of the receiver to the per- 
turbation of the sender can be very large. Close to the 
synchronization transition it diverges and even deep in- 
side the region of synchronization, where the transverse 
Lyapunov exponents are negative and large, huge excur- 
sions away from the synchronization manifold occur re- 
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FIG. 22: Panel (a) shows second moment xij (blue curves) 
and cross correlation C\j (red curves) for uniformly dis- 
tributed random noise, panel (b) shows bit error rate rij (blue 
curves) and cross correlation Cij (red curves) for binary ran- 
dom noise as a function of e. Combination of units ij as indi- 
cated in the legend. Network of four coupled Bernoulli maps 
with parameters a = 1.1, r = 100, a = 0.625 and M = 10" 8 . 



but the response can also be suppressed. 

Finally the linear response of a chain of three units and 
a network of four units has been investigated. We find 
that for a bi-directionally coupled chain of three units the 
perturbation is completely filtered out by the unit in the 
middle and both outer units can synchronize perfectly. 
Thus the second moment and the bit error rate becomes 
zero. 

The hub of a star network can receive simultaneously 
any number of secret messages without any error. 

In this work we restricted our investigations to chaotic 
maps which in some respects have different properties 
than chaotic flows. But with respect to complete syn- 
chronization, maps and flows are very similar and many 
of the obtained results are also observed in numerical 
simulations of chaotic differential equations. Hence we 
believe that our findings will contribute to a general un- 
derstanding of linear response of synchronized chaotic 
systems. 



Appendix A: Analytic Results For The Bit Error 
Rate 



Outside the synchronization region the bit error rate 
trivially equals 1/2. Inside the synchronization region in 
general, the bit error rate has to be determined by means 
of computer simulations. However, for some special cases 
it can be calculated analytically. 



1. Logistic map, uni-directional setup, r = 0, e = 1, 

K = 



suiting in a power law behavior and diverging moments 
of the distribution of deviations between the sending and 
receiving unit. Mathematically, this is a consequence of 
multiplicative and additive noise appearing in the equa- 
tions of linear response. 

The bit error rate of a transmitted binary message 
is used as another quantity to investigate the linear re- 
sponse. The bit error rate is given by an integral over 
the distribution of deviations between the trajectories of 
sender and receiver. For the unsynchronized system the 
bit error rate is at its maximum of 50%. Directly at the 
synchronization transition it decreases to smaller values 
and shows a complex nonmonotonic behavior inside the 
region of synchronization which cannot be related to the 
properties of transverse Lyapunov exponents. For special 
cases we could calculate the bit error rate analytically. 
Relating it to an iterated function system we found a frac- 
tal distribution of deviations yielding a devil's staircase 
for the bit error rate as a function of model parameters. 

The linear response to a periodic perturbation was in- 
vestigated as well. It shows resonances due to the delayed 
self-feedback of the sending unit. Depending on param- 
eters and frequency, those resonances can be very large 



In the case of the uni-directional setup with logistic 
maps and no time delay the bit error rate can be calcu- 
lated for the point e = 1. The dynamics is given by 



Xt+i = f(xt) 

Vt+i = f(x t -i + m t -i) ■ 



from which follows that 



H+l 



ft m t 



(Al) 



(A2) 



From the facts that (1) f' t and m t are uncorrelated, (2) 
the probability distribution of /' is symmetric about /' = 
and (3) m t — ±M, follows that d has basically the same 
probability distribution as /'. Therefore, the bit error 
rate can easily be calculated, see Eqs. (24 1 and (15): 



r = l(l ~ J P(d) ddj = Ml - I p(f)df j 



(A3) 



arcsec(4) 



0.4196 . 



This is in agreement with the numerical simulations. 
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2. Logistic map, bi-directional setup, r = 0, e = 4, 

K = 

Similarly, in the bi-directional case the bit error rate 
can be calculated for e = h. The dynamics is given by 



Vt+x = ^f{Vt) 



from which follows that 



1 



(A4) 



f(x t + m t ) 



1 



dt+x = ^/tW 



Comparing this with Eq. \A2\ leads to 




(A5) 



(A6) 



which is also in agreement with the numerical simula- 
tions. 



3. Bernoulli map, r = 



In Fig. 



11 



(a) one can discover a staircase structure for 
the bit error rate. For the uni-directional setup this is 
true for e > I, while for the bi-directional setup this is 



valid for 



< e < 



For these regions the bit error rate 



can be calculated analytically. If one takes a closer look 
at the staircases, Fig. [23] it becomes apparent that they 
have infinitely many steps, i.e. they are a kind of devil's 
staircase. 



a. Uni-directional coupling 

This staircase structure should be explained here for 
the case of uni-directional coupling. The equation for 
the distance can be written in the following way, see also 
Eq. (Ill: 



d. = 



4 



l(l-e)d t ^-\eM 
|(l-e)dt_i + |eM 



(A7) 



The two equations represent the two different bits. If d t 
is plotted versus dt-i, then d^ and d^~ are two paral- 
lel straight lines, see Fig. [24] The values of this itera- 
tion d t (d t -i) generate the distribution p(d) from which, 
in principle, the bit error rate can be calculated. The 
dashed boxes (- 



in Fig. 

dt is bounded to (due to t 
e < | the two maps d. 



24 represent the interval which 



re attracting fixed points). For 
and df have a certain overlap 
in their co-domain, see Fig. 24(a)| (the co-domains are 



indicated by gray stripes). This fact makes the distri- 
bution p(d) complicated. However, for e > | the two 



maps have no overlap, see Fig. 24(c) and the distribu- 
tion is manageable analytically. As a result of the gap 
in the co-domain (indicated by a zigzag pattern), a gap 
in the domain emerges in the next time step. The latter 
gap produces two further gaps in the co-domain which 
become gaps in the domain in the next time step. The 
result of this iterative process is that the distribution 
p{d) has a fractal support. The first and largest gaps are 
shown in Fig. [25] The first gap is called G. The gaps 
produced by G are G~ and G + . The gaps coming from 
G + are called G h and G ++ ; the gaps coming from G~ 
are called G and G H and so on. 

Now we want to calculate the exact position of the 
gaps. The fixed points of d~ and d + are called d~ and 
df. One can easily calculate that 



dt 



3eM 
3e- 1 

3eM 
3e- 1 



and 



From Fig. 24(c) it can be seen that 

G = ]d-(d+),d+(d-)[ 

3e(3e-2)M 3e(3e-2)M 



3e- 1 



3e- 1 



(A8) 
(A9) 

(A10) 



which is about ]-0.19M, 0.19M[ for e = 0.7, see Fig. [25 
The gap G + is generated by applying d + to G, i.e. 



G+ = }d + (d-(dt)),d + (d+(d:))[ 
3e(5 - 12e + 9e 2 )M 3e(5 



12e + 9e 2 )M 



2(3e- 1) 



2(3e- 1) 



(All) 



which is about ]0.96M, 1.14M[ for e = 0.7, see Fig. [25 

Then the gap G for example, is generated by ap- 
plying dr to G + , i.e. 



G~ 



}d-(d+(d-(d+))),d-(d + (d+(d-)))[ (A12) 



and so on. 

Due to the constant and equal slope of d~ and d + , 
and due to the fact that there is no overlap between the 
co-domains of d~ and d + , the relative frequency of all 
distances d which occur (i.e. which are not inside a gap) 
are equal. This means that Fig.[25]can also be seen as the 
corresponding histogram; all bars have the same height. 

The bit error rate is related to the integral from — M to 



+M over the distribution of the distances, see Eq. (24 1. 
From Fig. [25] it can be seen that for e = 0.7 this integral 
exactly equals |; thus, the bit error rate equals |> which 
is in agreement with Fig. [23] 

If e is changed, then the positions of the gaps are 
changed, too. As long as the gap G + contains the value 
+M (= as long as the gap G~ contains the value —M), 
the integral yields \ and the bit error rate is \. This 
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FIG. 23: Bit error rate r for uni-directionally coupled Bernoulli maps with r = and k — 0. Zooming in the staircase structure 
reveals more and more steps. 




(a)e = 0.6; d and d+ have an overlap (b)e = | (c)e = 0.7; d and d + have no overlap 

in their co-domains. in their co-domains. 



FIG. 24: The iteration dt(dt-i) = df(dt-i) of the distances for different coupling parameters e. The dashed boxes represent 
the interval which d is bounded to due to the fixed points. Additionally, the bisecting line dt(rft-i) = dt-i is plotted. The gray 
stripes show the co-domains of d~ and d + . 



explains the plateau AB in Fig. 
Eq. (All | we can calculate the 



plateau: 



2:-S 



With the aid of 
exact position of this 



A: 



B: 



3 



0.6, 

+ V2) 



0.804738 . 



(A13) 
(A14) 



Similarly, one gets the point C of Fig. [23] The bit error 
rate becomes | when the integral starts to be 0. This is 
when the gap G is as large as (or larger than) the interval 
[-M,M]. Considering Eq. ( |A10[ ) yields: 

C: e = -(3 + VE) « 0.872678. (A15) 



The plateau DE, which has the value |, can be calcu- 
lated considering the gap . One gets: 

D: e fa 0.837266 (A16) 
E: e w 0.866386 (A17) 

All other plateaus can be calculated with the aid of 
smaller gaps. 



b. Bi-directional coupling 

The calculations for the case of bi-directional coupling 
are very similar to the ones for the uni-directional case. 
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FIG. 25: Gaps in the domain/co-domain/distribution of d 
and d + for different recursion depths k. e = 0.7 . 



Here, only few results should be shown. 
For 



0.3 = \ < e < ^ = 0.5 (A18) 



the bit error rate is 0. 
For 



0.5749 w * + << | = 0.6 (A19) 
6 i/6 3 



the bit error rate is \. 
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